Phenomenological model for the normal state ARPES line shapes of high temperature 

superconductors 
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An adequate description of the single particle spectral function observed in high temperature 
superconductors in the normal state is an important goal, yet to be realized. Here, we show that 
with minimal phenomenological extension of the so-called extremely correlated fermi liquid (ECFL) 
model, it is possible to establish a phenomenological model, which is capable of describing ARPES 
data for different materials (Bi2Sr2CaCu208+(5 and La2-xSra:Cu04), with exceptional fidelity in 
fitting the data as function of momentum as well as energy. 
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In the sudden approximation theory [l| of the angle 
resolved photo-electron spectroscopy (ARPES), photo- 
electron counts, /(/c,co'), recorded as a function of mo- 
mentum (k) and energy (oo)^ are given by 

I{k,L0) = \Mif\^f{c0)A{k,L0) (1) 

where Mif is the dipole matrix element for the photo- 
excitation, f{uj) is the Fermi-Dirac distribution func- 
tion, and A{k^ u) is the single particle spectral function, 
A{k^uj) = ;^ImG(^, cj), where G is the single particle 
Green's function^. 

As the single particle Green's function in the normal 
state contains vital information on the nature of excita- 
tions relevant to the high temperature ( "high Tc" ) super- 
conductivity, its characterization by ARPES has been a 
major line of research. Various approaches towards get- 
ting at this information have been attempted: a phe- 
nomenological approach based on a simple scaling be- 
havior of the electron self energy Q , an asymptotic solu- 
tion to the Gutzwiller projected ground state of the t-J 
Hamiltonian [5] , application of a non- Fermi liquid theory 
[sl for low dimensions, and a newly proposed solution to 
the t-J Hamiltonian [7]. 

For an experimental "cut," i.e. an experimental data 
set taken along a hue of k values, I{k^uo) is a func- 
tion defined on a two dimensional domain. This multi- 
dimensionality makes analyzing I(k^uj) a non- trivial 
task. While attempts [8] have been made to analyze the 
I{k^uj) image (e.g., see Figure [3ja)) as a whole, the cur- 
rent theoretical understanding of line shapes in terms of 
A{k^ uj) depends on the analysis of selected energy distri- 
bution curves (EDCs; EDC means /(/c, uj) as a function of 
(J for a fixed value of k) [5, [H, 0, 0] or selected momentum 
distribution curves (MDCs; MDC means I{k,uj = ujq) 
with a constant ujq and k along a line) [ol, [lo| . 

Currently, there is no consensus on a theoretical model 
that can suitably describe ARPES data of high Tc mate- 
rials. A model that can describe the normal state data, 
both EDCs and MDCs, obtained in different experimen- 
tal conditions and for different materials with the same 



intrinsic parameters would be a good candidate. Here, 
we propose a new phenomenological model that satisfies 
just such a requirement. This model is built upon the re- 
cent success of the so-called extremely correlated Fermi 
liquid (ECFL) model R. This model retains all features 
of this past work, but with crucial extension that makes 
the model suitable for describing MDCs as well. 

The result is a phenomenological model in which an in- 
trinsic inelastic background ("caparison factor") due to 
extreme electron correlations 0, IHI continues to play a 
critical role in fitting EDCs while the effect of the capari- 
son factor on MDCs is significantly reduced, suggesting 
an interesting dichotomy between the behaviors of the 
EDC and the MDC, somewhat reminiscent of a similar 
dichotomy observed in low dimensional systems 0, [HI . 

A phenomenological study of this kind seems to be 
helpful, also in light of the on-going development of 
the ECFL theory [13, The theoretical formalism 

of ECFL initiated by Shastry [il|Jl3,] is quite involved, 
and, while a numerical solution valid for hole dop- 
ing X ^ 0.3 is now available, more time seems necessary 
to extend these promising results to near-optimal dop- 
ing. Thus, a phenomenological model based on the main 
feature of the theory, the caparison factor, may be of con- 
siderable value at this stage. In this theory [13], the ca- 
parison factor is an uj- dependent adaptive spectral weight 
that encodes two key pieces of physics: the Gutzwiller 
projection that reduces the spectral weight at high uj and 
the invariance of the Fermi surface volume at low uj. 

In our previous work [7| , it has been demonstrated that 
the normal state EDCs for optimally doped cuprates for 
two different compounds, or for different experimental 
conditions (low photon energy or high photon energy), 
can be explained using an ECFL line shape model, all 
with one set of intrinsic parameters. We will refer to that 
model as the "simplified ECFL (sECFL)" model [IF 




in relation to the fuller theory in development ^ ^ 

While the EDC analysis, the focus of Ref. 0, has strong 
merits [7|, ll6| , a natural subsequent question is whether 
MDCs can be described as well within the same theory. 
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FIG. 1. Line shape fits of EDCs for Bi2212 using (a) 
sECFL, (b) pECFL with momentum dependence (MD), and 
(c) pECFL with momentum independence (MI). Data and 
model parameters are identical with those in Ref. 7, except 
for additional parameters used for (b), ai = 1.2 and a2 = 2.5 
(cf. Eq.©, and a slightly different ry value (0.17 0.18 eV). 



In the sECFL model [7], the Green's function was 
taken to be of the form 



G{k,uj) 



uj - e(k) - $(cj) 



(2) 



1 - - 



= is the total spectral weight 



where Qr, 

per k in the t-J model, n is the number of electrons per 
unit cell, and x = 1 — n is the hole doping per unit cell, 
^{(jj) is an ordinary Fermi liquid self energy (Eq. 2 and 
footnote 6 of Ref. 7) containing two intrinsic parameters, 
ZpL and cjo, and one extrinsic parameter, rj. Aq is the 
energy scale that is uniquely determined by n^ZpL^^o 
through the global particle sum rule [l7| . 

The above Green's function can be rewritten as 




In UJ- e(k) - ^{uj) 
UJ — e{k) \ 



(3) 
(4) 



where Cn{k^uj) is the "caparison factor" 0, Ell and 7^ cx 
Aq is an energy scale parameter, to be discussed further 
below. As all symbols in Eq. [3] other than ^{uj) are real. 



A{k,uj) = Cn{k,uj) ApLik.uj) 



(5) 



where Afl is the spectral function for the "auxiliary 
Fermi liquid (AFL)" Green's function [18|, Afl = 
I ImGFL = ^ Im [a; - e{k) - ^{uj)]-\ 

We will show how the function Cn can be modified 
phenomenologically so that the main sECFL idea, sum- 
marized concisely in Eqs.|4]and[5l can be fully preserved 



while the range of applicability of the model can be ex- 
tended in a critical manner. We shall refer to the modi- 
fied model as pECFL, with p for "phenomenological." 

We begin our discussion by showing line shape fit re- 
sults first, and will then explain our new model. 

Figure 1 shows ARPES line shape fits for the normal 
state data for the optimally doped Bi2Sr2CaCu208+(5 
(Bi2212) sample along the "nodal direction," (0,0) 
(7r,7r). The data are identical with the data discussed 
in Figure 4(b) of our previous work (7|. Panel a shows 
the sECFL fit essentially [lo] identical with the fit in that 
work, while panels b,c show pECFL fits. With pECFL, 
the fit quality remains equally excellent for /c ~ /c^, while, 
for k away from kF , it has improved noticeably, albeit by 
small amount (e.g., see the EDO at top or bottom). 

Figure 2 shows ARPES line shape fits for MDCs of the 
same data set. Panel a shows clearly that sECFL has 
difficulty fitting the data even for a; = 0, i.e. the Fermi 
energy {Ef)- In panel b, where the pECFL(MD) model 
is used, the fit improves visibly. However, the quality of 
the fit of panel c, where the pECFL(MI) model is used, 
is clearly the best. 

We note that for the data fit in Figures 1 and 2 a small 
"elastic background line shape" (0.5 times the raw line 
shape for k ^ /ci?) had been subtracted, as explained 
in Ref. 7. Then, the data are fit fully with the spec- 
tral weight given by Eq. [T] without any additional back- 
ground intensity at all. That the model pECFL(MI) is 
able to describe both EDCs and MDCs in this manner, 
with completely identical parameter values, including the 
T] value, demonstrates the high fidelity of the model. In 
contrast, it is impossible to achieve optimal fits of MDCs 
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FIG. 2. Line shape fits of MDCs for Bi2212 using (a) sECFL 
and (b) pECFL(MD) and (c) pECFL(MI). Fit parameters are 
the same as those used for Fig. [Tl except for the 77 values for 
(a,b), as noted here. 



and EDCs with the same rj parameter within sECFL or 
pECFL(MD): as shown at top of panels a,b, the fit qual- 
ity degrades significantly if the same 77 value as that for 
EDCs is used. These conclusions remain the same for the 
fit of the 91 K data [T^fioj. 

Now, we explain our new model, and its two regimes 
("MD" and "MI"). 

From Eq. O it is required that Cnik^uo) > for any 
value of k and cj, or equivalent ly 7^ > u — e{k). For 
7n > 0, it can be seen that this condition is satisfied 
sufficiently close to the Fermi energy {uj = 0) and the 
Fermi momentum {€{k) = 0). On the other hand, for a 
negative value of e{k), uj—e{k) becomes increasingly more 
positive where important spectral features are expected: 
as the peak position is given by ~ Ze{k) with the mass 
renormalization ratio < Z < 1, we see that uj — £(k) 
grows as ~ (1 — Z)\e{k)\^ as eik) decreases from 0. 

In the sECFL model, 7^ is given by 4(5^Ao/n^, which 
we shall define as 7no- For our fit parameters [7], 
Ao = 0.12 eV, n = 0.85, we get jno = 0.38 eV. Since 
Zfl = 1/3 [7J and thus Z = ZflQu = 0.19, we see that 
Cn > requires that e{k) ^ —0.47 eV. For the data fit in 
FigurelU the smallest eik) value is —0.4 eV for /c = /ci, as 
shown in Figure [3fa). Thus, the simple truncation pro- 
cedure, (Cn max(Cn,0)), put in place for the sECFL 
in Ref. M to assure the positivity of the spectral function, 
was justified. 

How might the model be modified so that it is suit- 
able for describing MDCs at low energy as well, while it 
preserves, or improves, the quality of EDC fits? As we 
showed in the above figures, a dramatic improvement of 
MDC fits and a substantial improvement of EDC fits can 
be achieved at the same time. While the improvement 
of the EDC fit quality is arguably only modest for the 
fits shown in Figure [TJ it is clear from the above discus- 
sion that the improvement will be dramatic outside the 
momentum range shown in Figure [H Figure [Sfb) shows 
that this is indeed the case. 

Our new pECFL model is based on two generic obser- 
vations: (i) from a more fundamental point of view [21], 
going beyond sECFL, 7^ can be dependent on both £{k) 
and cj, and (ii) MDCs observed for high Tc cuprates tend 
to be fit well with Lorentzian line shapes in general (but 
see Figure H]). 

With these observations, the following can be noted. 
Considered as a function of uo — sik)^ Cn needs modi- 
fication only in the limit oi uo — e{k) ^ 00. If 7n is 
taken to be a function of uo and eik)^ this can be accom- 
plished by requiring that 1/7^ ^Oascj — £(^) 00. We 
keep 7n = 7n,o for the opposite limit uo — eik) —00. 
This minimal modification of the model has the addi- 
tional virtue of establishing the correct limiting behavior 
G{k^uj) Qn/^ for UJ 00, as well as for uj —00, as 
required by the sum rule fiij per k^ within the t-J model. 
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FIG. 3. (a) The ARPES data for Bi2212, used for previous 
figures. Rectangle E (M) marks the range of data fit in Figure 
[T]©. Circle symbols show e{k) values 0] used in the fits. The 
color scale represents ARPES counts that increase from green, 
blue (half max), white, to red (max), (b) The EDC fits for the 
k value corresponding to line V of panel a. (c) The caparison 
factors {Cn{k,uj)). In pECFL(MI), we take Cn{k = /cf,cj) 
of sECFL, or pECFL(MD), and apply it to any e{k). (d,e) 
The weights determined from MDC fits (d) and EDC fits (e). 
Rectangular shades mark the C(;-range (d) and the /c-range (e) 
corresponding to fits in Figures I1I2[ respectively. 



A simple way to implement the modification is to use a 
sigmoidal reduction of 1/7^ d^s uj — eik) ^00: 



7nO • 



1 + exp 



UJ - eik) - ai7no 



a27n0 



(6) 



The model that incorporates this 7^ function is named 
"pECFL(MD)," where MD means "with momentum de- 
pendence," as Cn{k,uj) is dependent [23| on k through 
e{k). 

The role of the two new parameters, ai and a2, is to 
provide "soft landing" to the caparison factor. With good 
values of ai and a2, the caparison factor never dips be- 
low zero, and approaches zero only asymptotically, from 
above, diS uj — e{k) 00. The key is that the soft landing 
must be engaged early (small ai), but not too early, and 
quickly (small a2), but not too abruptly. After some cal- 
culus, one sees that the success, i.e. the non-negativity of 
Cn for any value of uj — e{k), is guaranteed if 



ai < l + a2(l-loga2) ai,max(tt2) 



(7) 
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This condition was used as a constraint during the hue 
shape fit using the pECFL(MD) model. In practice, we 
find that only very narrow ranges of possible ai and a2 
values are allowed: both must be of order 1. 

Another regime for pECFL is motivated by the obser- 
vation (ii) above. Lorentzian-like MDCs can be taken 
to mean a momentum independent self energy [24]. This 
motivates us to consider simply dropping the s{k) term 
in Eq. HI and obtain a momentum-independent version 
of sECFL, which we cah "pECFL(MI)." However, since 
we are concerned only with ARPES fits here {uj ^ 0), 
pECFL(MI) can also be thought of as a momentum- 
independent version of pECFL(MD), whose caparison 
factor for e{k) = and ^ is practically identical 
with that of sECFL for the same parameters, as shown 
in Figure [3fc). 

Lastly, the quality of line shape fits can be examined 
by inspecting the weight parameter that determines the 
overall scale: Wm{^) = \M\'^f{uj) for MDC fits, shown in 
Figure [3fd), and We{k) = |Mp for EDC fits, shown in 
Figure (He). In this figure, these functions are shown for 
ranges of uj and k going beyond the ranges for Figures 
1 and 2. As the energy dependence of |Mp is expected 
to be weak for the narrow range of cu used in the fit, we 
expect that the shape of Wm{^) be determined mostly by 
f{oj) alone, which we show for comparison in Figure [3fd). 
Also, we expect that We{k) show only modest variation 
in the momentum range of the fit [H, |25| . It is clear that 
these expectations are met increasingly better as we go 
from sECFL, to pECFL(MD), and to pECFL(MI). For 
example, for sECFL, Wm{^) deviates strongly from f{oj) 
for UJ < —50 meV, and We{uj) increases by more than an 
order of magnitude as decreases. 

Therefore, for the Bi2212 data examined, pECFL(MI) 
emerges as the best model. However, that this may not 
be a universal behavior is indicated in Figure HI Here, 
the EDC data and the MDC data for La2-xSr^Cu04 
samples, as we found available in the literature 0,0, 
are fit. For the pECFL(MD), ai = 1.6 and a2 = 2.0 are 
used, while all other parameters are the same as those for 
Bi2212, except the different doping parameter n = 0.937 
necessary for panels d-f. For panels a-c, ujq = 0.42 eV is 
used, the value found [7] to fit the data to very high \uj\ 
value, compared to ujq = 0.5 eV, equally good for fitting 
data up to \uj\ ~ 0.25 eV, just as for Bi2212. 

Panels a-c of Figure 4 show that the EDC at kp 
is described well by all models [28], as we found for 
Bi2212. The MDC fits {uj = 0; panels d-f) however give 
an edge to pECFL(MD), compared to both sECFL and 
pECFL(MI): it gives half as small x^- Key observation 
is that the MDC dit uj = is notably asymmetric. 

In conclusion, in this Letter, we propose a phenomeno- 
logical ARPES line shape model based on the ECFL ap- 
proach. The essential characteristic of the model remains 
the caparison factor [7], which is modified in this Let- 
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FIG. 4. EDC fits to the LSCO data (a-c, optimally doped, 
n = 0.85), and MDC fits to the LSCO data (d-f, under doped, 
n = 0.937). The data for (a-c) are the same as those used in 
Figure 4(c-d) of Ref.0, while the MDC data in (d-f) are taken 
from Ref . 1271. 



ter to excellently describe both EDCs, characterized by 
strong asymmetry, and MDCs, seemingly more "conven- 
tional," much more so for Bi2212. While our model is 
not the first to fit both EDCs and MDCs (9| of high Tc 
cuprates, its demonstrated fidelity and range of appli- 
cability is now unprecedented. In addition, continuing 
first principle studies [l^ indicate that the interplay be- 
tween the extreme correlation and the one electron band 
structure (e.g., the t^t ratio) may be important - we 
hope to learn whether such interplay accounts for certain 
material dependence that we are beginning to recognize 
here. More experiments and line shape analysis, in co- 
ordination with theory development, may shed light on 
this and other remaining questions such as momentum 
dependence and doping dependence. 
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